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We describe a novel method which allows the treatment of high orders of replica-symmetry- 
breaking (RSB) at low temperatures as well as at T = directly, without a need for approximations 
or scaling assumptions. It yields the low temperature order function q{a, T) in the full range < 
a < oo and is complete in the sense that all observables can be calculated from it. The behavior of 
some observables and the finite RSB theory itself is analyzed as one approaches continuous RSB. 
The validity and applicability of the traditional continuous formulation is then scrutinized and a 
new continuous RSB formulation is proposed. 
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I. INTRODUCTION 

The ordered phase of spin-glass models [l| has gained 
much attention over the past three decades, but still 
there are many open questions even in models which 
were designed to possess simple 'mean-field- type' solu- 
tions. Paradigmatic for the difficulties encountered in 
the description of the ordered phase is the appearance of 
an ultrametric structure in the mean-field theory of spin 
glasses, known as Shcrrington-Kirkpatrick (SK) model 
[2] . The correct treatment of the SK model involves a hi- 
erarchical scheme as introduced by Parisi Q , which only 
recently has been proven Q to be exact in the limit of an 
infinite number of hierarchical steps of replica symmetry 
breaking (RSB). 

Though at temperatures right below the freezing tran- 
sition Tc, a small number k of RSB steps is not a bad 
approximation, at temperatures T <^ Tc the convergence 
of the results with respect to k becomes worse. Tradition- 
ally, the K = oo limit is formulated by a continuous the- 
ory where the functional free energy f[q{x)\ is maximized 
with respect to the Parisi order function q{x), < x < 1 
[U, [1]. This has later been expressed by a closed self- 
consistency form, involving a set of numerically solvable 
partial differential equations [1, Q. The treatment of 
RSB at finite temperatures is well established, mean- 
while. Near Tc, there are even analytical solutions avail- 
able. At zero temperature, however, the proper form of 
the theory is still under discussion. The standard differ- 
ential equations of the continuous theory and their initial 
conditions get singular in the zero temperature limit and 
the traditional formalism breaks down. 

In the literature, one finds two main directions of ap- 
proaching the zero temperature limit. Pankov proposed a 
beautifully simple scaling Ansatz directly at infinite RSB 
which becomes exact in a certain limit at zero tempera- 
ture 0|- This limit, however, is far from being the only 
important region of the full zero temperature solution as 
we shall show later. On the other hand, our group has 
developed an exact finite RSB approach, which involves 
some advanced numerics when trying to reach orders of 



RSB, high enough to obtain a confident extrapolation 
K — > oo. However, our approach allows for a numerically 
exact description of the full SK-modcl with all its diffi- 
culties and features at T = 0. The numerical results can 
be used to construct and test analytical approximations 
or simpler alternative theories which capture the essen- 
tial physics while being generalizable to more complex 
physical situations Q. 

The low temperature formulation, which is developed 
in the present work, directly identifies and analyzes the 
issues of the infinite RSB limit at zero temperature and 
resolves them by appropriate rcscaling of auxiliary quan- 
tities. It is based on the idea of rescaled Parisi block 
size parameters [1, H a, = Prrii with /? = T'^. By 
this transformation, the non-trivial structure of the zero- 
temperature order function of the SK-model can be re- 
solved. In the traditional formalism, the structure is con- 
centrated and hidden at the point a; = 0, where q{x) gets 
singular. Both, Pankov's and our approach, investigate 
this point, but while the scaling approach is incomplete 
in the sense that it e.g. only accounts for the short time 
observables or that it cannot describe the situation at fi- 
nite external fields correctly, our approach yields the full 
zero temperature solution. 

The paper is organized as follows. In section II, we give 
a short review of the issues encountered when performing 
the zero temperature limit and show the connections be- 
tween the two approaches mentioned above. In section 
III, we derive the non-trivial zero-temperature limit of 
the SK-model at finite RSB. In section IV, some results 
at low temperatures, obtained with the finite RSB for- 
malism are shown. In section V, we discuss the k ^ oo 
limit of our theory. 



II. THE T = LIMIT OF q{x) 

The Parisi order function q{x), x G [0, 1] which is 
the central quantity at finite temperatures becomes a 
constant function at zero temperature (limT_,o q{x) = 
1, a; > 0) with a singularity at a; = 0. Since typical 
observables as e.g. the internal energy u involve inte- 
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FIG. 1: (color-online): The true T = order function q{a,0) 
at 200 RSB (black) and the step approximation at 20 RSB. 
The lower (red) curve shows the result of Pankov-scaling, i.e. 
q{a, 0) = 1 — 7a~^. 



grals of the form it = /3 dx{l — q{x))'^, one is tempted 
to absorb the divergent factor /3 by a transformation 
X — *■ a = (3x. The non-trivial zero temperature struc- 
ture which was concentrated at a; = in the Parisi order 
function and created the singularity there is then 'blown 
up' and resolved by the new order function q{a). The 
transformation might remind the reader of PaT scaling 
(l3t . The idea here is different, though, as we don't seek 
for a universal form of the order function (which is known 
to be nonexistent) but for a non-singular zero tempera- 
ture theory. 

The transformation a; a is not limited to T = 
0, of course, so that we can formulate a general low- 
temperature theory which allows a smooth transition 
from a treatment directly at T = to finite tempera- 
tures. The objective of our formalism is to find the func- 
tion q{a, T) with a G [0, (3] from which one can obtain 
q{x,T) = q{Px,T) at finite temperatures. At k < 00, 
q{a, T) can be displayed as a step function with the num- 
ber of steps equal to k (see Fig. [T|). 

In this notation, the Pankov-scaling approach is equiv- 
alent to expanding q{a, 0) ~ 1 — 70"^ near a = 00 and 
calculating the coefficient 7. The knowledge of the a~^- 
term of this expansion is sufficient for short-time observ- 
ables as the entropy or the non-equilibrium susceptibility. 
The full solution, however, must respect the functional 
dependence of q{a,0) at all a G [0,00] which is needed 
to calculate e.g. the ground state energy or the equilib- 
rium susceptibility. Another point is, that the action of 
a finite external field is reflected by the small a domain 
of q{a,T) which is not resolved in Fig- [D the 

difference between the two solutions at low a is shown. 

In order to obtain a well defined zero temperature 
limit, we start with a formulation of the theory at fi- 
nite orders of RSB (k < 00), where q{a,T) is a function 
with K plateaus. The height of the ith plateau is given 
by the nunbcr qi G [0, 1] and the step positions are given 
by fli G [0,/?]. If one aims at a continuous theory for 
AC — > 00, one must first check that 5qi = qt — qi+i and 



Sai = fli — fli+i, i.e. the step heights and widths approach 
zero. For T > 0, those requirements are fuUfiUed and are 
the basis of the well known Parisi differential equations 
PiSH- At zero temperature and k — > 00, we still find 
Sqi — *■ which is not surprising, since qi is restricted 
to [0,1] independent of temperature. However, some step 
widths Sai grow indefinitely in this limit, so the assertion 
6ai — > which is needed to derive the Parisi differential 
equations fails at zero temperature. The non-vanishing 
Stti correspond to those which approach infinity for 
K — > 00 and so the continuous theory fails at a = 00, 
while for finite a, a differential equation formulation is 
valid. 

In section V of this paper, we show how the failure of a 
continuous theory at a = 00 can be overcome by changing 
the initial condition of a partial differential equation. 



III. RECURSION TECHNIQUE FOR 
ARBITRARILY HIGH RSB ORDERS k 

We investigate the standard SK-modcl Hamiltonian [l|, 
13] of iV Ising spins Si = ±1 with external field H 



HsK = ^ JijSiSj + Hy^ i 



(1) 



The quenched random coupling constants Jij are Gaus- 
sian distributed random variables with zero mean and 
variance iV~^[l3]- The model yields a freezing transition 
at Tc = 1- Since the coupfing constants' degrees of free- 
dom are quenched, the free energy per spin of the system 
is given by / = -T/N (log (exp{-PHsK)) s) j where {■) ^ 
refers to an average in the space of Ising-spins Si = ±1 
and {■) J refers to an average with respect to the cou- 
pling constants Jij. In order to bypass the average of 
the logarithm we avail ourselves of the standard replica- 
trick. After transformation to a single site model p^ . 
wc introduce Parisi RSB for the replica coupling matrix, 
perform the replica limit n ^ and the thermodynamic 
limit N —^ 00 and arrive at an expression for the free 
energy per spin (see also 



/ - - 2gi) + J ^ qnm^ - m,_i) - / (2) 



K+l 
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with K the order of RSB and qi the values of the elements 
of the Parisi-blocks with size rrii [3]. The non-trivial 
part / of the free energy is given by a k + 1-fold nested 
integral 
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where two Gaussian integral operators are defined as 
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f{hi) 
(4) 
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with Ag; = q., - qi+i, Ag^+i = Qk+i and 



GE 



G 



(5) 



The last Gaussian integral in this chain is centered at the 
external field ft.„+2 = H. 

As explained above, it is convenient at low temper- 
atures to transform to new variables = (3mi. Since 
^ _ _ fH+i_ definitions of the Gaussian integral 

operators / and / remain invariant. The explicit 
temperature dependence is located at the initial condi- 
tion, which has a well defined zero temperature limit in 
the a-formulation 



(2cosh(/3/ii)) 



(6) 



The utility of our new formulation for low temperatures 
lies in the non-singular limit of this inner integrand. It 
has been shown before @ that the first integral can be 
performed analytically in the T ^ limit. For finite 
temperatures, there is also an asymptotic regime hi 
±oo of the integrand where 

(2cosh(/3/ii))^"i e'^il'^il. (7) 

In the asymptotic regime, the Gaussian integrals can 
be performed analytically at each level of integration 



GE 



exp ( a; ( \hi+i \ + l-a^Aqt 



(8) 



while the deviation of the integrals from its asymptotic 
form at small h is incorporated in form of auxiliary func- 
tions expCj(/i) defined by 



GE rG 



(2cosh(/3/ii)) 



Tai 



exp ^a; ^ fliAgj + \h.,+i\ + expCj(/ij+i)| ^ 



(9) 



The functions expCj are well defined and free of singu- 
larities at each level of integration, for all temperatures 
including T = directly and at all k < oo. It is therefore 
the object of our choice for the finite RSB numerics. At 
zeroth level of integration the 'initial condition' 



expCo(/i) = Tlog 1 + exp 



2N 
T 



(10) 



has a smooth transition from finite temperatures to T = 
0; i.e. liniT^o cxpCg(/i) = 0. Furthermore, the expC 
obeys a h —i- —h symmetry at each level of integration - 
irrespective of temperature T or magnetic field H and is 
continuous at ft, = 0, though there is a kink at /i = with 
9/iexpCj(/i)|^^Q+ = —1, Vi. Therefore, it is sufficient to 
restrict ourselves to /i > and introduce a boundary 
condition at /i = in the continuous theory. 

At finite k, there are simple recursion relations between 
expCj(ft-) and cx.pC^_i{h) 



J 



expC,(/i) = — log 



exp ( a, ( \h'\ - \h\ + expC,_i (ft') 
I 



■^a^Aq^ 



(11) 



and in the end. the non-trivial part of the free energy can 
be expressed in terms of an integral of expC^ 

/ = (expCJft) + |ft|) + i ^ akAqk (12) 
•^''+1 k=i 

so the free energy per particle can be written as 



/ = -jxL + ^E«. [(g. - 1)^ - (gm - If] 



G 



(expC Jft) + |ft|) (13) 



K + l 



with the non-equilibrium susceptibility Xne — /3(1 — ^i). 
For the continuous RSB limit, one can fix the domain 
boundaries of g(a, T) by defining a^+i = and ao = 



f3 and the first line of becomes proportional to an 

integral da{l — q{a))'^ . 

At finite n, the free energy must be maximized with 
respect to {ai,qi}. This is done for T = (finite T) 
by a root search of the gradient of / in the 2k {2k + 1) 
dimensional parameter space jlTj . 

There are two main advantages of performing numerics 
in terms of expC. The first is that the non-trivial domain 
of this function (i.e. where it is non-zero, see Fig. [5]) is 
strongly restricted to the region of small ft and so the 
numerical effort reduces considerably. This is illustrated 
by Fig. [3] where the K-dependence of typical calculation 
times is shown. A single calculation implies the calcu- 
lation of the free energy and all its derivatives. Since 
the number of derivatives grows as 2k -I- 1, the full calcu- 
lation time grows stronger than linear. The calculation 
time normalized to the number of derivatives to calcu- 
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initial condition given in equation (jlOp . 
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FIG. 2: Non-trivial domain of expCf{h) of a typical 50-RSB 
calculation at T = and H — represented by the bars in 
h direction. The non-trivial domain is defined as the interval 
where expCj(/i) is numerically greater than zero. 
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FIG. 3: Calculation times at T = in dependence of k. 
The calculation time of the free energy with all derivatives 
(left part) goes with k^'^ and the calculation time of a single 
derivative (right part) goes with k^'^ 



late, however grows nearly linearly with k and therefore 
high orders of RSB may be obtained. The second advan- 
tage is, that expC(/i) ~ 0{1) and no loss of numerical 
precision due to large numbers happens. 

Figure [4] shows the functions expCj at all integration 
levels in case of zero temperature and for finite tempera- 
tures. The initial condition is represented by the lowest 
curve, which is exactly zero for T = and remains fi- 
nite for finite temperature corresponding to the non-zero 
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FIG. 4: (color-online): expC^{h) at even levels of integrations 
at T = and n = 50. The lower lines (red) correspond to 
small i and the upper lines (blue) to large i. The inset shows 
a T = 0.1 calculation with k = 30. 



IV. RESULTS 

In this section, we discuss some results of our exten- 
sive numerical computations according to the scheme de- 
scribed in section III. We have been able to perform cal- 
culations at extremely high orders of RSB (up to 200 
at T = and up to 53 at finite temperatures) with un- 
precedented accuracy which allow a deep insight into the 
subtleties that arise as k approaches its physical limit 
K = oo. 

In a previous publication we used 42 RSB data to ex- 
trapolate the zero temperature free energy to k = oo 
by fitting to the function = + ^|22^[g. 

The extrapolation is now confirmed according to 200 RSB 
high precision numerics with the new expC-formalism to 
/^^^ = -0.763 166 726 566 547, a = 4. It is difficult to 
estimate the error of the extrapolated free energy. The 
fluctuations of the free energy of the finite k calculations 
arc on an arbitrary small scale, due to the utilization of 
arbitrary precision arithmetics, so that the accuracy of 
/oo depends only on the fit procedure. To give the reader 
an idea of the convergence level at 200 RSB, let us state 
that /200 — /oo — 10~^^. With standard fit methods, an 
accuracy of 10""'^^ is obtained, which represents the min- 
imum accuracy of /oo . With some technical tricks which 
are beyond the scope of this paper, however, it seems 
that the accuracy might be up to 10~^^. In any case, 
this is by far the most precise ground state energy ever 
obtained for the SK-modcl and it provides a good test for 
all coming formalisms which work directly at k = 00 and 
T = 0. In comparison with the literature, we find that 
the value is consistent with the estimate of Parisi p^ . 
but is slightly outside of the confidence interval given in 
i. 

At zero temperature the free energy / equals the in- 
ternal energy u. At finite temperatures, they are still 
closely related. We have checked numerically, that the 
free energy can be expanded in a Taylor series with reg- 
ular exponents at T = 

/(T) = /o - soT + hT^ + hT^ + 0{T^) (14) 

where sq = — Xne/4 is the zero temperature entropy, 
which can also be expressed in terms of the non- 
equilibrium, or single- valley susceptibility Xne- The in- 
ternal energy can be expressed by the same coefficients 
u{T) = f -T% = fo- /2T2 - 2hT^ + 0(T4), and 
also the temperature dependence of the entropy can be 
written as s(T) = sq - 2/2T - 3/3r2 + 0(T''). From our 
results we find /2 ^ and /a —0.24, in agreement 
with the literature 0, 0] . 

The zero temperature entropy is directly related to the 
non-equilibrium susceptibility which can be written at 
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finite RSB as Ul 



aiiqf 



li+i) 



(15) 



One finds, that Xne oc (k + kq) with 7 = 1.666664 ± 
5 • 10^^. For an exponent, this numerical accuracy is 
sufficient to claim 7 = 5/3. As a result, both, the non- 
equilibrium susceptibility and the entropy vanish with 
irregular exponents 5/3 and 10/3 for k ^ 00 at zero tem- 
perature. The implications of those irregular exponents 
will be discussed elsewhere [T2| . 

As we have explained above, the natural formulation 
of the order function at low temperatures is in terms of 
a function q{a,T) since it resolves the structure at the 
singular point x = of the original Parisi order function 
T) at T = 0. In order to better understand the crit- 
ical properties of the zero temperature order function, 
we have to discuss the finite RSB step approximation to 
q{a,T). In Fig. [5] the logarithm of the Parameters are 
plotted as a function of log k. One can clearly see that 
for large the spacing between successive a-points does 
not vanish, while for moderate a^, a continuum emerges. 
Also at small a;, a spacing appears on a loga-scale. This 
small-log a spacing, however, does not imply a spacing 
on an a scale, where the differential equations of a con- 
tinuous theory appear. At finite temperatures, the dis- 
creteness at large a disappears due to the restrictedness 
of Oj to the interval [0,/3], while the small a discreteness 
remains. A finite external field would destroy also the 
small-a discreteness. The notion of non-zero plateau- 
widths on a log scale is directly related to the discrete 
spectra of the rj-levels, which will be discussed more thor- 
oughly elsewhere [l^ . 

At T > 0, the identification of a break point x is im- 
portant. One finds, that 



lim Tax. 



(16) 



From our computations at finite temperatures, we can 
extract a confident value for x at temperatures down to 
T = 0.015. Below this temperature, calculations at k > 
50 are needed. There is, however, no reason to expect a 
large deviation of x{T = 0.015) from its true extrapola- 
tion to T = 0. We find x{T ^ 0.015} = 0.54684±2 • 10-^ 
in consistence with the literature 

UM- At zero temper- 
ature, X is ill defined because q{x,0) = 1, Vx > 0. As a 
consequence, one does not find x directly in the T = 0- 
theory. In fact, there is a subtle non-commutativity of 
the r — > and k —^ 00 limits for quantities like the 
break-point [T3 |. 

At finite k and T = 0, no m,; remains finite, all go to 
zero. From scaling arguments, however, one can see that 
there are finite rrii even at T = in the k = 00 limit. 
These are exactly the rrii = Tai corresponding to the 
discrete anoints (see Fig. [S|) and so the finite rrii are also 
discrete [20| . In the continuous formulation of the infinite 
RSB limit, an initial condition of a partial differential 
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FIG. 5: log ai as a function of log n for n = 1, .., 200 at T = 
and for k = 1, .., 50 for T = 0.03. 



equation is given at mi [2l| . There is, however, the region 
of discrete rrii where a continuous theory is invalid, so at 
T = the initial condition of the continuous theory is 
disconnected from the validity domain of the differential 
equation. This is why the traditional theory fails for 
T = 0. 

Two other important and closely related parameters 
are the T^-coefficient of the T-expansion of the Edwards- 
Anderson order parameter qEA{T) = 1 - aT^ -I- 0{T^) 
and the a~'^-coefficient of the 1/a-expansion of g(a,0) = 
1 - 7a~^ -I- 0{a~^). At finite orders of RSB, there is 
also a linear term in the T-expansion of qEA with the 
coefficient equal to Xne, but this coefficient vanishes for 
K ^ 00. From our numerical data, we can extract a = 
1.59411 ± 0.00002 and 7 = 0.4108 ± 0.0001. 

In the following paragraph, we analyze our results from 
the viewpoint of PaT scaling [l3| . In the original formu- 
lation, PaT scaling has been used to evaluate a (approx- 
imately) universal scaling function q{x,T) ~ f{x/T) for 
X < X. This scaling is known to not hold exactly. If 
it held exactly, then f{x/T) = q{x/T,0) at all tempera- 
tures. In Fig. [HI we compare q{a, 0) to the result of the 
PaT hypothesis, which was obtained according to the de- 
scription in [l3| . 

In order to discuss the correction to PaT scaling near 
T = 0, we write 



q{a,T) 



q{a,0) 
q{a,0) 



q{a,T) 
q{a,T) 



for 
for 



a < a 
a > a 



(17) 



with a — j3x and g(a, T) is the correction to PaT scaling 
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V. oo-RSB LIMIT 

As K — > oo, the number of parameters {qi,ai) goes to 
infinity and a smooth order function q(a,T), defined on 
[0,/3] rises for T > 0. In this hmit, the set of functions 
expCj become one continuous function of the variables a 
and h defined as 



expC{ai,h) = Um cxpC,;(/i) 

K — ^OO 

and the free energy for H = can be written as 



(25) 



FIG. 6: (color-online): Comparison of the universal scaling 
function obtained by the PaT-hypothesis (black, dashed) and 
the zero temperature 200 RSB order function q{a, 0) (red, 
solid). The inset shows the difference of the two curves. 



near T — 0. From the identity 1 = J^^ da{l — q{a, T)) one 
can derive the correction to the PaT-scaling break point 
at zero temperature 



1 1 d 



(18) 



with a the quadratic temperature coefficient of the 
Edwards- Anderson (EA) order parameter qEA{T) ~ 1 — 
aT^ . If we assume that q{a,T) can be expanded in a 
Taylor-series near {a,T) ~ (oo,0), i.e. 



7(a,T) =^T*q,(a), for T « 1 



and 



9i(a) = ^ b{o- ^ , for a > 1 



(19) 



(20) 



one can derive relations between the lowest coefficients 
of the q{a, T) expansion, the quadratic temperature co- 
efficient a of the EA order parameter, the break point 
and the a"^ coefficient 7 of the expansion of g(a,0) at 
a = 00. 



6? = 6} = 



/ daqi{a) = ^(1-25;) -6^. 
Jo X 



(21) 
(22) 

(23) 



The first relation states that qi (a) must go to zero faster 
than as a ^ 00. The parameters 7, a, i are very well 
known from the literature and have been obtained from 
our numerical data, too. From the above relation, one 
can thus extract = -0.22035 ± 0.00012 and write the 
correction to PaT-scaling near {a,T) = (00, 0) as 

q{a, T) = q{a, 0) - 0.22 + 0(T, a~^f. (24) 

Relation can be used as a check for qi (a) extracted 
from numerics. 



1 

f = -- da{q{a) - if - expC(a = 0, /i = 0) (26) 
4 Jo 

where expC(a,ft.) is obtained by solving the partial dif- 
ferential equation 



daexpC = - 



q{a) 



(9^cxpC + 2a9/iexpC + a (9/iexpC)' 



with boundary and initial conditions 

5/iexpC(a, /i)|^^g+ = —1; cxpC(a,oo) = 
expC(/3, h) = Tlog (1 + exp {~2(3h)) 



(27) 



(28) 



for h>0. 

At finite temperatures, where the spacing between suc- 
cessive ai and qi go to zero as k — > 00, the expC is 
well-behaved and the differential equation can be solved 
numerically. In this case equations (j27l28p are merely 
a reformulation of the Parisi theory [l], |3[, convenient 
for investigations at low temperatures. At exactly zero 
temperature, however, two problems arise which force 
us again to reformulate. Firstly the discreteness in the 
large a regime of q{a) formally invalidate the differential 
equation at a = 00, where the initial condition is given. 
Secondly the second derivative of the initial condition has 
a divergence at = as T — *■ 0. 

The second issue can be resolved by a further rescaling 
y = {a + l)h and introduction of the function g{a,y) ~ 
(a -|- l)expC(a, y/{a + 1)). The initial condition at (3 for 
expC translates to g{P,y) ~ log(l -|- e~^^) = go{y)- The 
left part of Fig. [7]shows the a and y-dependence of g{a, y) 
(for visualization purposes, the a-axis has been rescaled 
to the interval [0,1] by the introduction of the variable 
C = a/(l + a)). It shows that for finite a (i.e. C, < 1), 
g{a, y) varies smoothly as a function of a and the descrip- 
tion in terms of a partial differential equation is valid. In 
the large a limit, however, a singularity appears. To un- 
derstand this singularity and the solution of this issue, we 
go back again to the finite RSB formulation. The inves- 
tigation of the 200 RSB calculations shows that at large 
a, the difference between g{ai, y) and g{ai^i, y) does not 
vanish with k ^ 00 as it does for moderate a <^ ai (2^ 
(see right side of Fig. [7]). This means, that the recursion 
relation (jlip docs not approach a differential equation at 
a = 00. Instead, one can see that the recursion drives 
the function to a limiting function g^o (y) from which the 



7 





0.2 0.4 0.6 0.8 1 



^ a+1 

FIG. 8: (color-online): g{C,0) with ( — a/ {I + a) for up to 
200 orders of RSB. The dots with the different colors refer to 
low (red) and high (blue) RSB orders. The solid line (red) is 
obtained by the oo RSB formalism proposed in the text. 



continuous part of g{a, y) starts. In the right part of fig- 
ure [7] the first 10 results of the recursion relation for a 
K = 200 calculation are shown together with 50 (y) ^md 
.900(2/), which is obtained by numerically solving an ordi- 
nary differential equation, as explained below. 

To complete our discussion, we now approach the point 
a — 00 from below. In the finite a regime, the behavior 
of g{a,y) is governed by the partial differential equation 
(the dot refers to a derivative with respect to a, while the 
prime means a y-differentiation) 

g = -% + !) [{a + l)g" + 2ag' + aig'f] + 

2 a+1 

(29) 

In order to investigate this equation at zero temperature 
in the limit a — > 00, where the initial condition is given, 
we expand the order function q(a,0) = 1 — 70"^ and 
equation itself near a = 00. To first order in we 



find g = ^F[g] with 



F[g]^g-yg'-j{g" + 2g' + {gr)- (30) 



Obviously, an initial condition g(oo,y) of (|29p with 
F[g[oo,y)] 7^ would lead to a logarithmic singularity 
of g{a,y) at a = 00. The only non-singular initial con- 
dition is therefore the solution to the ordinary differen- 
tial equation F[goo{y)\ = 0. Indeed, the solution g{y) of 
this differential equation seems to be the limiting func- 
tion goo (y) of the recursion starting from go (y) discussed 
above. In Fig. [7l the function goo{y) is the numerical 
solution of F[goo] = 0. In some sense, the partial dif- 
ferential equation governing the function g(a,y) yields 
its own initial condition - it is the only initial condition 
which makes sense. 

To further illustrate and confirm this line of reasoning 
and to better understand the transition k ^ 00 we shall 
restrict our discussion of g{a, y) to y = as representa- 
tive for the a-dependence of g. In Fig. [8]we plot g(ai, 0) 
at T = for K < 00 varying from 10 to 200. Again, 
one can see the discreteness at a = cx). For demonstra- 
tion purposes, a numerical solution of (|29p is plotted as 
the drawn through line. To obtain this solution at large 
a, equation ((29|l has been expanded up to order at 
a = 00. The solution of this expansion taken at a = 8 
(C ~ 0.89) is then used as an initial condition for the 
full partial differential equation (|29|) at a = 8. By thor- 
oughly looking at the line, one can see a small error near 
this junction point. With more effort like higher order 
expansions at a = 00 or advanced numerical methods for 
partial differential equation (e.g. pseudo-spectral meth- 
ods), the quality of the full continuous RSB solution at 
zero temperature can be strongly improved, but this is 
beyond the scope of this work. 
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VI. CONCLUSION 

We have developed an RSB technique which allows cal- 
culations at extremely high orders of RSB near and at 
T = 0. We have indeed performed calculations at T = 
for up to 200 RSB and for finite temperatures up to 53 
RSB. We extracted the dependence of various observables 
on the order of RSB and on temperature and obtained, 
to our best knowledge, the by far most precise numerical 
value for the ground state energy of the SK model. The 
connection to PaT-scaling has been discussed and a cor- 
rection to the PaT-scaling assumption q{x,T) = f{x/T) 



near zero temperature and a = x/T = oo has also been 
discussed. Furthermore, we have proposed an Ansatz for 
a full treatment of the zero temperature limit of the SK- 
model directly at infinite RSB - in analogy to the con- 
tinuous RSB formalism at finite temperatures. It would 
be usefuU to derive a closed set of self-consistency equa- 
tions for q{a) in the sense of Sommers and Dupont Q 
and solve them numerically in order to obtain the zero 
temperature order function directly in its physical limit. 

We gratefully acknowledge useful discussions with 
David Sherrington, Kay Wicsc and Markus Miiller. 
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